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We investigate the pairing physics of a three-component spin-orbit coupled Fermi gas in two 
spatial dimensions. The three atomic hyperfine states of the system are coupled by the recently 
realized synthetic spin-orbit coupling (SOC), which mixes different hyperfine states into helicity 
branches in a momentum-dependent manner. As a consequence, the interplay of spin-orbit coupling 
and the hyperfine-state dependent interactions leads to the emergence of Fulde-Ferrell (FF) pairing 
states with finite center-of-mass momenta even in the absence of the Fermi-surface asymmetry that is 
usually mandatory to stabilize an SOC-induced FF state. We show that, for different combinations of 
spin-dependent interactions, the ground state of the system can either be the conventional Bardeen- 
Cooper-Schrieffer pairing state with zero center-of-mass momentum or be the FF pairing states. 
Of particular interest here is the existence of a three-component FF pairing state in which every 
two out of the three components form FF pairing. We map out the phase diagram of the system 
and characterize the properties of the three-component FF state, such as the order parameters, 
the gapless contours and the momentum distributions. Based on these results, we discuss possible 
experimental detection schemes for the interesting pairing states in the system. 

PACS numbers: 03.75.Ss, 03.75.Lm, 05.30.Fk 


I. INTRODUCTION 

Ever since its experimental realization, synthetic spin- 
orbit coupling (SOC) in ultracold atomic gases has at¬ 
tracted much attention BHT]. Over the past few years, 
SOC-induced exotic phases and phase transitions have 
been extensively studied in both the Bose and the Fermi 
gases [HHISl- Central to the effects of SOC in these 
systems is the nonperturbative modification of the single¬ 
particle dispersion spectra. For Fermi gases in particular, 
the SOC-modified single-particle dispersion can induce 
exotic few-body states [niiiH] as well as highly nontriv¬ 
ial many-body correlations [inillO]- The SOC-induced 
Fulde-Ferrell (FF) pairing state is an interesting example 
where the interplay of SOC and the Zeeman-field-induced 
Fermi surface asymmetry stabilizes unconventional pair¬ 
ing superfluids with finite center-of-mass momenta m- 

Recently, it has been shown that an alternative FF 
pairing mechanism exists in a Fermi-Fermi mixture where 
a two-component noninteracting Fermi gas, while dressed 
by SOC, interacts spin selectively with a third fermionic 
species [20]. As the hyperfine-spin distribution in the 


* qinfang@phy.ccnu.edu.cn 
t wzhangl@ruc.edu.cn 
t wyiz@ustc.edu.cn 


SOC-induced helicity branches is asymmetric in momen¬ 
tum space, the spin-selective pairing interaction leads to 
pairing states with finite center-of-mass momenta in both 
the two-body and the many-body sectors. This novel 
pairing mechanism is thus fundamentally different from 
the majority of existing proposals of SOC-induced FF 
states since it is originated from spin-selective interac¬ 
tion and does not require a Fermi-surface asymmetry. 

In this paper, we extend this exotic spin-selective- 
interaction-induced FF mechanism to a three-component 
Fermi gas where all three hyperfine states are cou¬ 
pled by SOC HU From the single-particle disper¬ 
sion, we show that in the three SOC-induced helicity 
branches, two exhibit asymmetric hyperfine-spin distri¬ 
butions. Based on this asymmetry, we discuss several 
different configurations of the spin-selective interactions, 
under which the ground state of the system can either 
be a conventional Bardeen-Cooper-Schrieffer (BCS) state 
with zero center-of-mass momentum or be FF states orig¬ 
inating from the interplay of SOC and pairing interac¬ 
tion. In particular, an interesting three-component FF 
state can be identified in which every two out of the 
three hyperfine components form FF pairing with a com¬ 
mon center-of-mass momentum. We investigate the sta¬ 
bility and phase transitions of the three-component FF 
states by mapping out the mean-field phase diagram. De¬ 
pending on the excitation gap of quasiparticles, the FF 
states in the present system can be further categorized 
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FIG. 1. (Color online) Level diagram of the Raman processes 
generating the synthetic spin-orbit conpling. The Raman 
laser beams are counterpropagating along with freqnen- 
cies uji and uji -\- AtJi (i = 1,2). The two-photon detuning is 
5, and e is the quadratic Zeeman shift m- 


as gapless or fully gapped. We characterize the gapless 
FF states by studying the gapless contours in momen¬ 
tum space and discuss their relation with features in the 
number distribution, which should facilitate experimen¬ 
tal detection. As the various configurations of interaction 
can in principle be experimentally implemented via the 
Feshbach resonance technique, our paper not only reveals 
the generality of the FF pairing mechanism induced by 
SOC and spin-selective interactions, but also has interest¬ 
ing implications for future experiments on SOC-induced 
exotic superfluidity. 


The paper is organized as follows: In Sec.|^ we calcu¬ 
late the single-particle dispersion under SOC and analyze 
the inherent asymmetry in the hyperfine-state distribu¬ 
tion. In Sec. m we present the system as well as the 
mean-held formalism. We then consider different conhg- 
urations of spin-dependent interactions and investigate 
their impact on the pairing superfluidity in Sec. |IV| In 
Sec. |Vj we discuss the novel three-component FF state. 


The main results of this paper are summarized in Sec. VI 


II. SINGLE-PARTICLE DISPERSION 


We consider a two dimensional three-component Fermi 
gas where the three atomic hyperhne states are labeled 
as I ± 1) and |0). The hyperhne states are coupled by Ra¬ 
man lasers as illustrated in Fig. Under an appropriate 
rotating frame, the single-particle Hamiltonian can be 



FIG. 2. (Golor online) (a) Single-particle energy spectra along 
the fey = 0 axis in momentnm space for a homogeneous gas. 
The red dashed curve denotes Ei, the blue dashed-dotted 
curve is E 2 , and the black solid one represents E 3 . (b)-(d) The 
momentum distributions of the hyperhne states in different 
helicity branches. The dashed curves denote |aip, the solid 
curves are |feip, and the dashed-dotted ones represent |cip. 
Here, i = 1 — 3 for (b)-(d), respectively. The parameters are 
he = 0.44E,., hS = 0, and h = 2.425E,. [H]. 


written in momentum space as [HllIIlEs] 


Hoik) = 

{k+2kj.e^)'^ 


2m 


h 

0 


— he 


0 

h 

{k-2kr-e^)^ 

2m 


+ hS 


( 1 ) 


where S is the two-photon detuning of the Raman pro¬ 
cess, e accounts for a small quadratic Zeeman shift, and 
h = hil.ji/2 denotes the strength of the effective Zee- 
man field, which is proportional to the Rabi frequency 
of the Raman process We take the recoil energy 

Er = /{2m) and the corresponding wave vector k^ 

as the units of energy and wave vector, respectively, with 
m as the atomic mass. In the rest of the paper, we will 
only consider the case where 5 = 0. The Hamiltonian 
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Ho(fc) can be diagonalized as 


Ho{k) 

a* (A:) 

= E,{k) 

ai{k) 

biik) 

h{k) 


. Ci{k) _ 


. Ciit) . 


where the eigenvalues Ei{k) (i = 1 — 3) are the single¬ 
particle dispersion spectra for the helicity branches, and 
= {ai,bi,Ci)’^ are the corresponding eigenvectors. 
Importantly, the coefficients ai{k), bi{k), and Ci(fc) are 
related to the weight of hyperfine states in the corre¬ 
sponding helicity branches. 

In Fig. we plot the single-particle dispersion spectra 
as well as the momentum distributions of the hyperfine 
states in the helicity branches along the ky = 0 axis. Ap¬ 
parently, the hyperfine-state superpositions in the helic¬ 
ity branches are momentum dependent [Figs. [^b)-[^d)]. 
A critical observation is that although the momentum 
distribution of the hyperfine state |0) is symmetric with 
respect to = 0 for all three helicity branches, such 
a symmetry is absent for the cases of states |1) and 
I — 1). Instead, the momentum distributions of |1) and 
I — 1) are symmetric with respect to each other such that 
= \ci{—kj;)\'^. The presence or absence of this 
inversion symmetry has crucial effects on pairing physics. 
In fact, in the weak-coupling limit, pairing tends to oc¬ 
cur between two fermions residing on their corresponding 
Fermi surfaces in the absence of SOC. With SOC mixing 
up the hyperfine spins into the helicity branches, should 
we turn on a small attractive interaction between hyper¬ 
fine states |1) and | — 1), the pairing state would have 
zero center-of-mass momentum. On the other hand, if 
we turn on a small attractive interaction between states 
|1) and |0) (or | —1) and |0)), the pairing state would have 
a finite center-of-mass momentum. This is similar to the 
FF pairing mechanism in Ref. [20] where the interplay of 
SOC and spin-selective interactions plays the key role. 


III. MEAN-FIELD TREATMENT OF PAIRING 
STATES 


With the understanding of the single-particle disper¬ 
sion, we now study the many-body pairing physics using 
the standard mean-field formalism. Now the complete 
Hamiltonian can be written as 


H = J2 <i>Hk)[Ho{k) - fi]^k) + Hi, 

k 


( 3 ) 


where is the creation operator 

for the three hyperfine states, ^ is an annihilation oper¬ 
ator of a Fermi atom with wave vector k and pseudospin 
i = 1,0, —1 describing three atomic hyperfine states, ct 


k.i 


is the corresponding creation operator, and fi is the chem¬ 
ical potential. Here, the most general form of the inter¬ 
action Hamiltonian H/ is 




ct - C,-/ nC/I 1 
k,l Q-k,0 .0 Q-k',1 


+ 


Q,k,k' 
ffl-1 
V 


E 


ct ct. - Cc, IT/ 1 

k,l Q-k,-l k',-1 Q-k',1 


Q,k,k' 


, 50,-1 t t 

+ V- E Q-k-l^k'-l^Q-k',0- 

Q,k,k' 


Here, Q is the center-of-mass momentum of the pairing 
fermions, V is the quantization volume in two dimen¬ 
sions, and gij {i,j = ±1,0) is the bare interaction rate 
between different hyperfine states, which can be renor¬ 
malized following the standard procedure in two dimen¬ 
sions [24H26] . 


1 _ 1 V ^ 

gi,j V ^ hP-k'^/ra + Eb^ij ' 


( 5 ) 


Here, E^^ij is the two-body bound state energy between 
states |i) and \j). The two-body bound state energy can 
be related to the corresponding three-dimensional s-wave 
scattering length, which can be tuned experimentally via 
the Feshbach resonance technique. 

Under the mean-field approximation, the interaction 
term becomes 


k 

+E -1 ^3.1-i^±i) 

k 

+ E (^Q.O,-l4,o'^Efc-1 -i^±o) 

^|Aq,i,_iP 

5i,o 5i.-i 5o,-i 

where the superfluid order parameter is taken as 


= (bj = l,0,-l). (7) 


Therefore, in the hyperfine-spin basis 'l'^(fc) = 

- .crnict; - ,ct i,ct. - Y'. the effective 
V fc.i’ fc.o’ Q_fc,o’ fe.-i’ Q-fc.-U ’ 

mean-field Hamiltonian can be rewritten as 


2 E'I'! 3 (*)^feEW 

k 


E ( n^Q-k + 


- V 




5i,o 




Q,0,-ll 
5i,-i 5o,-i 


( 8 ) 


The matrix takes the form of 


( 4 ) 
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FIG. 3. (Color online) Contour plots of the thermodynamic potential landscape in the Ai,o-Qj; plane for (a) /i = lOF^r 
and Ebi = O.SF^r; (b) = 5Er and Ebi = 0.3Er’, (c) ^ = 2Er and Ebi = O.SEr. The global minima are located at (a) 
(Ai,o/F^r ~ 2.3410, Q:,/kr ~ -2.0184); (b) (Ai,o/F^r ~ 1.5917, Q^/K ~ -2.0347); (c) (Ai,o/£^r ~ 0.5970, Q^/kr ~ -2.0666). 
Other parameters are chosen as he = 0.2Er, M = 0, and h = Er. 
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with /(2m) — fi. 

The zero-temperature thermodynamic potential can 
then be derived 
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T-)-0 




fc,j = l,2,3,4,5,6 

^3 




- V 


|A, 


Q.i.ol 


|A, 


Q.l-ll 


|A, 
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5i.o 


51-1 


5o,-i 


( 10 ) 


where 0(a;) is the Heaviside step function, Tr denotes 
the trace over both the momentum and the spin degrees 
of freedom, T represents the temperature, and fcs is the 
Boltzmann constant. 

In the following, we first consider the simple case where 
inter-atomic interactions are limited between two spin 
species. Due to the symmetry of the setup at 5 = 0, we 
only need to study two different cases among all three 
possible combinations of interactions: the case with in¬ 
teraction between |1) and |0) and the case with interac¬ 
tion between |1) and | — 1). 


IV. PAIRING STATES WITH INTERACTION 
EXISTING BETWEEN TWO SPIN SPECIES 


We first study the case where interaction only presents 
between states |1) and |0) so that A^^ , = A^^ o -i “ 

0. We also define A^ ^ ^ = Ai_o and E},i = Eh,i,o to 
simplify notations. Under these conditions, we calculate 
the thermodynamic potential in Eq. (101 for any given 
Q by numerically diagonalizing the matrix Eq. (|^. For 
the parameters that we have studied, the local minimum 
in the thermodynamic potential always occurs with Q = 
Qx^x-: where is the unit vector along the direction of 
the SOC. In Figs. |^a)-[^c), we show the typical contour 
plots of the thermodynamic potential on the plane of 
^ 1 , 0 -Qx for a given chemical potential /i and binding 
energy E{,i. From these figures, we can see clearly that 
in general there exists only one local minimum, which 
corresponds to the ground state of the system. This is 
consistent with our previous analysis that under SOC- 
induced asymmetric hyperfine spin distribution and spin- 
selective interaction between |I) and |0), the pairing state 
acquires a nonzero center-of-mass momentum, i.e., it is 
an FF state. 


These observations allow us to minimize the thermo¬ 
dynamic potential with respect to the order parameter 
Ai 0 and the center-of-mass momentum Q = Qj. to find 
the ground state of the system. In Fig. |^a), we map out 
the phase diagram on the /r-Fti plane for he = 0.2Er, 
hS = 0, and h = Er. Apparently, a continuous phase 
boundary exists between the FF state, which is charac¬ 
terized by a finite Ai^o and a finite Qx, and the normal 






















5 



FIG. 4. (Color online) (a) Phase diagram on the ^-Ebi plane 
for he = Q.2Er, hS = 0, and h — Er- The blue solid 
curve represents the phase boundary of the continuous phase 
from the FF state to the normal state, and it marks the 
|Ai,o| = 10“^i5r threshold. The blue dashed curve denotes 
the phase boundary between the fully gapped FF state and 
the gapless nodal FF state, (b) Typical ground state super¬ 
fluid order parameter Ai,o and center momentum of the 
pairing fermions versus jj, for he = 0.2Er, hS = 0, h = Er 
and Ebi = Er. Typical contours of gapless points in momen¬ 
tum space for the gapless nodal FF state with (c) /r = 20Er, 
Ebi = 2Er (nFF); (d) /i = 5Er, Ebi = 2Er (nFF). 
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state (N), which is characterized by a vanishing order 
parameter. Furthermore, judging from the minimum ex¬ 
citation gap, both a fully gapped FF state (gFF) and a 
gapless nodal FF (nFF) state exist on the phase diagram, 
which are separated by a continuous phase boundary. In 
Fig.ib), we show how the ground-state order parameter 
Ai 0 and the center-of-mass momentum Q^. evolve with 
the chemical potential fi with fixed h = Er and E^i = Er¬ 
in Figs. Be) and B^)> we also show the typical gapless 
contours of the nFF state in momentum space. 

For the case where the only interaction in the system 
is between states |1) and | — 1), we may set ^ ^ = 
0-1 = Fy diagonalizing the effective Hamiltonian 
and minimizing the thermodynamic potential as we have 
performed previously, we find that the ground state of the 
system is either a BCS pairing state with Q = 0 or a nor¬ 
mal state with a vanishing order parameter. This is also 
consistent with our previous analysis in the weakly inter¬ 
acting limit that attractive interactions between states 
|1) and I — 1) would lead to BCS pairing due to the sym¬ 
metry in the momentum distribution of these states in 
the helicity branches. 


FIG. 5. (Color online) (a) Phase diagram on the fi-Ebi plane 
for he = 0.2Er, hS — 0, h — Er, and Eb 2 = Er- The red solid 
curve is the first-order boundary, whereas the green (light 
gray) solid curve represents the phase boundary of the con¬ 
tinuous phase from the FF state to the normal state, and the 
blue dashed curves denote the phase boundaries between the 
fully gapped FF state and the gapless nodal FF state. The 
black dashed line denotes Ebi = Er, and it is the asymptote 
of the first-order boundary, (b) Typical ground state super¬ 
fluid order parameters Ai,o, Ai,_i and center momentum 
of the pairing fermions versus /r for he = 0.21?^, hS = 0, 
h = Er, Eb 2 = Er, and Ebi = 1.15i?r. The red solid curve is 
Ai,o, the blue dashed curve represents Ai,_i, and the black 
dashed-dotted curve denotes Qx- 


V. PAIRING STATES WITH MULTIPLE 
INTERACTIONS AMONG VARIOUS SPIN 
COMBINATIONS 

In this section, we consider the case where attractive 
interactions are not only present between states |1) and 
|0), but also exist between states |1) and | —1). Hence, we 
can set A,^ g = 0, A^^^ = Ai _i, and Eb 2 = 

for convenience. 

As before, we numerically minimize the thermody¬ 
namic potential to look for the ground state of the sys¬ 
tem. In Fig. Ba), we show the typical phase diagram 
on the ^-Ebi plane for he = 0.2Er, h6 = 0, h = Er, and 
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FIG. 6. (Color online) Contours of gapless points in momentum space for different novel gapless nodal FF states with (a) 
^ = 20Er, Eh\ = 2Er (nFFl); (b) p = 5iJr, Eb\_ ~ 2Er (nFFl); (c) /r = 20Er, Ebi ~ Q.2Er (nFF2); (d) /i = 5i?r, Ebi = Q.2Er 
(nFF2). (e)-(h) Quasiparticle (quasihole) dispersion spectra along the ky = 0 axis for the novel gapless nodal FF states that 
correspond to panels (a)-(d), respectively. Other parameters are chosen as he = 0.2Er, hS = 0, h = Er, and Eb 2 = Er- 




k /k 




k /k 


FIG. 7. (Color online) Number density distributions in momentum space along the fcy = 0 axis for different novel gapless nodal 
FF states with (a) /r = 20Er, Ebi = 2Er (nFFl); (b) = 5Er, Ebi = 2,Er (nFFl); (c) fj, = 20Er, Ebi = 0.2Er (nFF2); (d) 

= 5Er, Ebi = 0.2Er (nFF2). The black dashed-dotted curve is tik^ilrir, the blue dashed curve represents nkfi/^r, and the 
red solid curve denotes nfc,_i/nr. Other parameters are chosen as he — 0.2Er, hS = 0, h = Er, and Eb 2 = Er- The unit of 
density is dehned through Ur = krl{2iT). 


Eb 2 = Er- It can be clearly seen from Fig.j^a) that there 
are four distinct FF phases on the phase diagram, includ¬ 
ing two nodal (or gapless) (labeled as nFFl and nFF2) 
and two fully gapped (gFFl and gFF2) ones. Whereas 
the gapped FF states are separated from the nodal FF 
states by continuous phase boundaries [blue dashed lines 
in Fig.j^a)], different nodal FF states or different gapped 
FF states are separated by a first-order phase boundary 
[red solid line in Fig. [^a)[. In this phase diagram, we 
also identify a normal state (iV) by numerically setting a 
small threshold of |A| = 10“^i5r- 

By adjusting the binding energy Ebi, the system can 
be tuned across different phases over a wide range of /r. 
In Fig.[5];b), we show how the ground-state order param¬ 
eters Ai.o, Ai__i and the center-of-mass momentum 
evolve with the chemical potential /i with fixed h — Ey, 
Eb 2 = Er, and Ebi = For this typical parame¬ 

ter set, the system can successively go through multiple 
phase transitions by increasing /i. In the local-density 
approximation where the effect of a global trapping po¬ 
tential is taken into account by the spatial variation in 


chemical potential, this is the order of phases that one 
would observe starting from a trap edge to its center. 

To further characterize the properties of these FF 
states, we demonstrate in Fig. the typical gapless con¬ 
tours and dispersion spectra in momentum space. For the 
gapless contours, we find that there are two closed gap¬ 
less rings in nFFl, and the two rings are both symmetric 
about the ky = 0 axis [see Figs. a) andj^b)). By de¬ 
creasing the chemical potential /i, the two rings become 
smaller in size, and gradually separate from each other. 
The gapless contours of nFF2 are shown in Figs. I^c) and 
l^d). Different from the two gapless rings in nFFl, there 
is only one gapless ring in nFF2 with given parameters. 
This gapless ring is symmetric with respect to the origin. 
As the chemical potential decreases, the gapless ring be¬ 
comes smaller. We also show in Figs. [§;e)-[§[h) the quasi¬ 
particle and quasihole dispersion spectra of nFFl and 
nFF2 along the ky = 0 axis. These are consistent with 
the corresponding gapless contours in Figs. ||;a)-[|[d) . In 
principle, one may probe these dispersion spectra ex¬ 
perimentally using momentum-resolved radio-frequency 
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spectroscopy. 

The interesting features of the gapless contours and the 
dispersion spectra also leave signatures in the particle- 
number distribution in momentum space, which may be 
probed more directly via the time-of-flight measurement. 
We show in Fig. the number density distributions in 
momentum space along the ky = 0 axis for different nodal 
FF states with the same parameters as in Fig. For 
the nFFl case with Ei,i = 2Er, it can be clearly seen 
from Figs.j^a) and[^b) that abrupt changes are present 
in the momentum-space density profiles, particularly for 
the I — 1) state. Comparing to the corresponding results 
of gapless contours, one can find that the discontinuous 
features are consistent with the right half of the gapless 
contours. For the nFF2 phase as shown in Figs. [^c) and 
[^d), similar discontinuities can be found in the momen¬ 
tum distribution of the |0) state, which are consistent 
with the corresponding structure of gapless contours. 

VI. SUMMARY 

We have studied the properties of the SOC-induced 
FF pairing states in a two-dimensional three-component 
Fermi gas with the recently realized synthetic SOC at 
zero temperature. The FF state here is the result of 
SOC-induced asymmetric momentum distribution of hy- 
perfine states and the spin-selective interaction, both of 
which are experimentally achievable. To illustrate this, 
we investigate in detail the impact of different combina¬ 
tions of spin-selective interactions on the properties of 
the pairing states of the system. Interestingly, the inter¬ 
play of SOC and spin-selective interaction can give rise 
to a novel three-component FF state in which every two 
of the three components form an FF pairing state with a 
common center-of-mass momentum. We study in detail 


the stability region of the FF states, the dispersion spec¬ 
tra of quasiparticle and quasihole excitations, the gap¬ 
less contours and the number distributions in momen¬ 
tum space, and discuss possible experimental detection 
schemes based on our results. As both the synthetic SOC 
and the spin-selective interactions have been realized ex¬ 
perimentally, our paper has interesting implications for 
future experiments on SOC-induced exotic superfluid¬ 
ity. In particular, we stress that even at temperatures 
above the superfluid transition temperature, this spin- 
selective pairing mechanism would assist the emergence 
of two-body bound states under appropriate interaction 
parameters. These two-body bound states should then 
acquire a finite center-of-mass momentum depending on 
the details of the interaction. As a consequence, the ex¬ 
otic pairing physics discussed above can be verified via a 
radio-frequency spectroscopy analysis of two-body bound 
states in the normal phase. 
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